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Abstract 



In view of the recent interest in reproducing holographically various properties of confor- 
mal fluids, we review the issue of vorticity in the context of AdS/ CFT. Three-dimensional 
fluids with vorticity require four-dimensional bulk geometries with either angular momen- 
tum or nut charge, whose boundary geometries fall into the Papapetrou-Randers class. The 
boundary fluids emerge in stationary non-dissipative kinematic configurations, which can 
be cyclonic or vortex flows, evolving in compact or non-compact supports. A rich network 
of Einstein's solutions arises, naturally connected with three-dimensional Bianchi spaces. 
We use Fefferman-Graham expansion to handle holographic data from the bulk and discuss 
the alternative for reversing the process and reconstruct the exact bulk geometries. 



Review to appear in the proceedings of the Corfu Summer Institute 2011: School and work- 
shops on elementary particle physics and gravity September 4-18, 2011, Corfu, Greece. 
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1 Introduction 

The relationship of fluid dynamics with general relativity goes back to the work of Damour 
PX Lately a different and perhaps more concrete shape of this relationship has been given 
by the so-called fluid/gravity correspondence (see for example the recent review 021) - Ac- 
cording to the latter, the gravitational degrees of freedom that reside in the boundary of an 
asymptotically AdSo spacetime describe the hydrodynamics of a relativistic fluid in D — 1 
dimensions. Consequently, the dynamical equations of the latter systems (e.g. Euler or 
Navier-Stokes) are encoded in the asymptotic behaviour of the bulk Einstein equations. 

The fluid /gravity correspondence framework appears to be capable of describing dif- 
ferent facets of relativistic fluids, such as superfluidity and dissipation. This appears to 
be a novel mechanism of emergence in physics whereby the low-energy effective degrees 
of freedom arise holographically in the boundary of a gravitational system. At a practical 
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level, fluid/ gravity correspondence currently occupies a large part of the AdS/ CMT corre- 
spondence as are collectively called the efforts to find new computational tools for strongly 
coupled condensed matter systems using holography (see e.g. the reviews HHSl)- 

Despite some important results in the study of holographic fluids, the issue of vortic- 
ity has been less well understood. This is important if one wants to extend the realm of 
AdS/CMT to interesting condensed matter systems such as rotating Bose or Fermi gases 
EEL turbulence or wave propagation in moving metamaterials (e.g. |8|). With such ex- 
tensions in mind, we review here our recent attempt to setup a holographic framework for 
the description of fluids with vorticity. Even without touching the thorny question of dis- 
sipation, i.e. assuming local equilibrium and non-dissipating kinematics, our studies reveal 
a remarkably rich structure as soon as vorticity is switched on. In particular, we note the 
intimate relationship of our neutral rotating holographic fluids with charged fluids in mag- 
netic fields, as well as with the problem of wave propagation in moving media. We believe 
that the latter observation can lead to the holographic description of analogue gravity sys- 
tems EHI21. 

In the present review we choose to devote most of our discussion to the gravitational 
side of the duality and we extensively discuss in Sec. |2] relativistic fluids, in Sec. [3] the gen- 
eral stationary Papapetrou-Randers and Zermelo geometries and in Sec. 0] various concrete 
examples of stationary geometries in 2 + 1 dimensions. In Sec. [5] we review the Fefferman- 
Graham construction of holographic fluids in 2 + 1 dimensions focusing on the Kerr-AdS4, 
the Taub-NUT-AdS4 and the hyperbolic NUT-AdS4 solutions. In Sec. |6]we give an overview 
of some new results, to be presented elsewhere 1151 , that aim to connect our approach with 
alternative descriptions of holographic fluids. Section [7| contains our conclusions. 

2 Relativistic fluids 

In this section, we recall the salient features of relativistic fluid dynamics (see e.g. lfl4llT5l ). 
This includes aspects of vector-field congruences and properties of the energy-momentum 
tensor. We work in arbitrary spacetime dimension D. 

2.1 Vector-field congruences 

We consider a manifold endowed with a spacetime metric of the generic form 

ds 2 = g llv dxf'dx v = Vab e a e b . (2.1) 

We will use a,b,c,... = 0, 1, . . . , D — 1 for transverse Lorentz indices along with a, /3, 7 = 
1, . . . , D — 1. Coordinate indices will be denoted ]l,v,p,.,. for spacetime x = {i, x) and 
i,j,k,. . . for spatial x directions. The dual of the orthonormal coframe e" is the frame e a , 
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which satisfies e a (e b ) = 5?. To define parallel transport we take the Levi-Civita connection 
coefficients Y? defined via the spin-connection one-form co a b as 

de» + cv\Ae b = 0, <V b = T&r, V h h = T c ab e c . (2.2) 

Consider now an arbitrary timelike vector field u = u a e a , normalized as rj ab u a u h = — 1, 
later identified with the fluid velocity. Its integral curves define a congruence which is char- 
acterized by its acceleration, shear, expansion and vorticity: 



wit 



V a u b = -u a a b + — — - ©h ab + cr ab + co ab (2.3) 



a a = u h V h u a , ® = V a u a , (2.4) 

= \K c K d ' {VcU d + V d u c )--^—^h ab h cA V c u d (2.5) 

1 

= V (fl M 6) +a (a u h) - — — -h ab V c u c , (2.6) 



D-l 

i 

"ab = ( VcMrf ~ Vd " c ) = V [a tl b] + a \a u b\ (2-7) 

The latter allows to define the vorticity form as 

2co = co ab e a A e h = du + u A a . (2.8) 
These tensors satisfy several simple identities: 

u a a a = 0, u a a ab = 0, u a co ab = 0, u a V b u a = 0, h c a V b u c = V b u a . (2.9) 

Killing vector fields, satisfying V( fl ^u = 0, are congruences with remarkable properties. 
We quote two of them, the proof of which is straightforward: 

• A Killing vector field has vanishing expansion. 

• A constant-norm§ Killing vector field is furthermore geodesic and shearless. It can 
only carry vorticity. 

The timelike vector field u can be used to decompose any tensor field on the manifold 
in transverse and longitudinal components with respect to itself. The decomposition is per- 



^ur conventions are: A (fl6) = 1/2 {A ab + A ha ) and A[ gfc ] = 1/2 (A ab - A ba ). 

2 This is not an empty statement since Killing vectors cannot be normalized at will. When their norm is 
constant, it can be consistently set to —1, or +1. 
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formed by introducing the longitudinal and transverse projectors: 



U\ = -u a u b , h\ = u a u b + 5 a b , (2.10) 

where h a y is also the induced metric on the surface orthogonal to u. The projectors satisfy 
the usual identities: 

U\U\ = U\, U a c h c b = Q, h a c h c b = h\, U\ = \, h\ = D-l. (2.11) 

For example, any rank-two symmetric tensor T ab can be decomposed in longitudinal, trans- 
verse and mixed components: 

T a b = eu a u b + S ab — u a q b — u b q a , (2.12) 
the non-longitudinal part being 

^ab = Sat — u a q b — u b q a . (2.13) 

We have defined 

e = u a u h T ab , S ab = h a % d T cd , q a = h a h T bc u c , (2.14) 

such that 

u a q a = 0, u a S ab = 0. (2.15) 

Finally 

u a T ab = q b - eu b . (2.16) 
2.2 The energy-momentum tensor 

The mere existence of a metric and a timelike vector-field congruence does not necessarily 
imply the presence of a relativistic fluid. If we wish to identify the timelike vector u with the 
velocity of a relativistic fluid, then we should require the presence of an additional symmet- 
ric rank-two tensor field - the energy-momentum tensor T ab whose projection along u is the 
(positive) energy density e of the fluid 

T ab u a u b = e, (2.17) 

as measured in the local proper frame. The latter concept deserves a comment. In non- 
relativistic fluids, the velocity field is unambiguously defined as the velocity of the mass 
flow of the fluid. In the relativistic case, it requires a more formal definition as energy and 
mass cannot be distinguished, and energy flows can be the result of dissipative phenomena 
or thermal conduction. One way to define the velocity, which amounts to defining a specific 
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local proper frame known as Landau frame, is to demand the absence of mixed terms in (I2.13I I: 



u a Z ab = 0. (2.18) 

Let us continue applying the decomposition (|2.12| l to the energy-momentum tensor. In- 
serting (12.131 1 in (I2.18L Eqs. (I2.15I I imply that q a vanishes in the Landau frame, where the 
energy-momentum is thus 

T ab = eu a u b + S ab . (2.19) 

The last piece S a b is the stress tensor, purely transverse. 

For a perfect fluid, all information is encapsulated in a further unique piece of data: the 
pressure p measured in the local proper frame. Hence, the stress tensor reads: 

S V f- = ph ab . (2.20) 

For a viscous fluid, the stress tensor contains friction terms: 

Sab = Phab + t a h ( 2 -21) 

where t ab is usually expressed as an expansion in the derivatives of the velocity field. At 
lowest order one finds 

t ab = -2 Wab - lh ab © (2.22) 

with n, £ the shear and bulk viscosities. In 2 + 1 dimensions there exists another term at 
this order, breaking the parity symmetry: Ca £ c M a uC(7 ^},y The coefficient £h is the rotational 
Hall viscosity. It characterizes a transport phenomenon similar to the Hall conductivity of 
charged fluids in magnetic fields. 

The dynamical equations for the fluid (Euler, Navier-Stokes, . . . ) are all encoded in the 
covariant conservation of the energy-momentum tensor 

V fl T flb = 0. (2.23) 

In the non-relativistic limit, Eq. (I2.23I I also delivers a matter-current conservation, which, for 
relativistic fluids, must be introduced separately as a consequence of charge conservation, if 
any. 

2.3 Effectively perfect fluids 

Relativistic fluids in the hydrodynamic regime are long wavelength approximations of finite- 
chemical-potential and finite-temperature states of certain (unknown) quantum field theo- 
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ries§l Quite generically all such fluids exhibit dissipative phenomena as they describe media 
with non-zero shear viscosity. However, all such fluids can be in special kinematic configura- 
tions where the effects of dissipation are ignorable^. In this case, their dynamics is captured 
by the perfect part of the stress tensor and the equations of motion reaco 

(e + p)0 + V l7 e = / s 

v r (2.24) 

(e + p)a — V±p = 

(V j_ = V + u V,-, stands for the covariant derivative along the directions normal to the veloc- 
ity field). Under this assumption, taking also into account the conformality (e = (D — l)p<x 
T D ), Eqs. (IZ241 lead to 

V i7 e = 

(2-25) 

a = ^F- 

The energy density is conserved along the fluid lines and in the absence of spatial pressure 
gradients (i.e. for energy and pressure constants in spacetime), the flow is geodesic. 

In several instances, the velocity of the fluid turns out to be a Killing vector field. Then, 
from the discussion of Sec. l2.1l several straightforward conclusions can be drawn: 

• The flow is geodesic, shearless and expansionless. 

• The internal energy density is conserved and the pressure is spatially homogeneous. 

• If the fluid is conformal then e = (D — l)p oc T D is constant in spacetime. 

Therefore, despite its viscosity, the kinematic state of the fluid can be steady and non-dissipa- 
tive. For this to happen, however, the existence of a constant-norm timelike Killing vector is 
required. In other words, the background geometry must itself be stationary^. In this case, 
the constant-norm timelike Killing vector congruence allows for the definition of a global 
time coordinate, with associated inertial frames. The latter are comoving with the fluid. All 
the examples that we will discuss in the following fall into this class. 



3 A special class of such fluids, actually the one that naturally arises in holography, are conformal fluids, 
i.e. those having vanishing energy-momentum trace: (D — l)p — e = (D — 1)£0. This equation is supposed to 
hold for any kinematic configuration, in particular when the fluid is at rest, where e = (D — l)p. The latter is 
therefore adopted as a thermodynamic equation of state valid always locally. When the fluid is not at rest, we 
conclude then that = 0, which must hold for any 0. In this scheme, the bulk viscosity for a conformal fluid 
is thus vanishing identically. Similar conclusions are reached for higher-order viscosity coefficients entering the 
traceful part of the energy-momentum tensor. 

4 A fluid can be stationary and altogether dissipate energy provided it is not isolated. These situations are 
better designated as forced steady states. On curved boundary backgrounds, the forcing task can be met by gravity 
through the boundary conditions. This was discussed in fl6l . As this feature does not appear in the backgrounds 
that will be analyzed in the forthcoming sections, we will not pursue this further. 

5 The interested reader can find more information on these specific issues in e.g. I17II18I . 

6 More generally, it can be shown that the velocity field m' 1 of a stationary fluid flow has to be proportional to 
a Killing vector field of the background geometry (171 . 
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3 Papapetrou-Randers stationary geometries 

Starting with appropriate time-independent bulk backgrounds^, conformal fluids appear 
holographically, evolving generally on stationary but not necessarily static boundary ge- 
ometries. Those fluids possess therefore non-dissipative dynamics inherited from the grav- 
itational environment and this dynamics contains in general vorticity. We will present here 
some basic properties of the boundary backgrounds arising in this context and explain how 
they affect the fluid dynamics. We postpone to Sec.|5]the actual holographic analysis relating 
some of these backgrounds to exact bulk Einstein spacetimes. 

3.1 General properties, geodesic congruences and Papapetrou-Randers frame 

Stationary metrics appearing in the holographic analysis we will be presenting later on are 
of the generic form 

ds 2 = B 2 (-(dt - bidx 1 ) 1 + a ij (x)dx i dx') , (3.1) 



where B, bj, a^ are x-dependent functions. These metrics were introduced by Papapetrou in 
||T9ll . They will be called hereafter Papapetrou-Randers because they are part of an interesting 
network of relationships involving the Randers form |20| , discussed in detail in H2TH and 
more recently used in [22U23I . 

For later convenience, we introduce , V and 7 such that 

"S* = 4. *-*h -> 2 = T^hkFf (3 - 2) 

The metric components read: 

goo = -B 2 , g 0l = B 2 bi, g lj = B 2 {a lj -b i b j ) / (3.3) 
and those of the inverse metric: 

? oo = L_ p o; = V_ ij = (3 4) 

& j2 B 2' & B 2 ' & B 2- y°-*> 



Finally 



B D y/a, (3.5) 



where a is the determinant of the symmetric matrix with entries fl,y. 

In the natural frame of the above coordinate system {d t , 3,}, any observer at rest has nor- 
malized velocity u = -^dt and dual form u = —B(dt — b) (we set b = bidx 1 ). The normalized 



7 Notice that in some instances time independence is not met in the bulk, but stationarity remains valid on the 
boundary as a consequence of appropriate boundary conditions 1 161 . These cases belong to the class of forced 
dissipative steady-states mentioned in footnote|4] 
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vector field u is not in general Killing - as opposed to d t . For this observer, the acceleration 
is thus non vanishing: 

a = Van = g ij di In B (dj + bjd t ) . (3.6) 

As already mentioned, the motion is inertial if and only if B is constant. It will be enough 
for our purposes to consider the case B — 1, and all subsequent formulas will assume this 
choice. We will furthermore introduce a frame 

So = d t , i H = Ej(bid t + di), E^E\ = 6{ (3.7) 

adapted to the geodesies at hand and its dual coframe (orthonormal as in Eq. (|2.1[)) 

e° = dt - b, t = E^dx 1 , E" f E^ = dij. (3.8) 

This will be referred to as the Papapetrou-Randers frame. 

The constant-norm Killing vector field il = Sq = d t (with u = —e° = —dt + b) defines 
a geodesic congruence (the orbits of all observers at rest in the Papapetrou-Randers frame). 
As was shown in Sec. 12.31 the latter has zero shear and expansion, but non-trivial vorticity (see 
Eqs. (27), dZ8l ): 

1 1 

co = -db => cv 0i = 0, cvij = - (dibj - djbi) . (3.9) 

The physical effect of vorticity is seen in the obstruction to the parallel transport of the spatial 
frame e a along the congruence: 

V lo g u = colfS^e y & V dt d i = to ij ai k {d k + b k d f ) (3.10) 

(tOij given in (|3.9|) are the spacetime components of the vorticity, while ca™ = EjEjcdij are its 
components in the Papapetrou-Randers frame). Embarked gyroscopes undergo a rotation. 

Papapetrou-Randers metrics do not exhibit ergoregions since§ goo = — 1. However, re- 
gions where hyperbolicity is broken (i.e. where constant-f surfaces become timelike) are not 
excluded. This happens whenever there exist regions where bibja') > 1. Indeed, in these 
regions, the spatial metric a^ — b[b- } possesses a negative eigenvalue, and constant-f surfaces 
are no longer spacelike. Therefore the extension of the physical domain accessible to the in- 
ertial observers moving along u = 3 f is limited to spacelike disks in which bjbja 1] < 1 holds. 
We will come back to this important issue in Secs. |3.3l and |5.2[ 

Before moving to the next topic, we would like to make a last remark. Following Eqs. 
(I2.25D , the shearless and expansionless geodesic congruence under consideration could de- 
scribe the fluid lines of a dissipationless stationary, conformal fluid, under the assumption 
that energy (and pressure) be conserved and constant all over space. As we will see in Sec.|5j 

8 Ergoregions would require a conformal factor in j3.H that could vanish and become negative. 
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this is exactly the dynamics that emerges through holography. 



3.2 Zermelo frame 

Trading the data (an, bj) for (hij, W) defined as 

hi = ^ L , V k h kl =5), (3.11) 



W = - 7 2 b\ Wi = huW = -h, (3.12) 



ij — y 2 , n n kj — Oj 

uW = 
J y 

the Papapetrou-Randers metric (I3.1D can be recast in the following form 

ds 2 = 7 2 [-d£ 2 + hij (dx* - W'df ) (dx> - W*'df) ] . (3.13) 

The latter is called Zermelo metric because it first appeared in the framework of the Zermelo 
problem M24H , yet another member of the relationship network mentioned above§. 

The Zermelo form of the metric suggests the following orthonormal coframe and its dual 
frame: 

z° = 7df / z K = L a ( (dx'- W'df), ^1^ = 5^, (3.14) 
z = \ (3t + W%) , z a = Ljdu = !%■ (3-15) 

We will call the latter the Zermelo frame. Its timelike vector field Zq defines a congruence of 
accelerated lines (Vz Zo 7^ 0). Thus, this frame in not inertial. It is instructive to compare the 
Papapetrou-Randers frame introduced previously (in d3.7| |, d3.8l l) with the Zermelo frame at 
hand. Being both orthonormal, they are related by a local Lorentz transformation, as one 
sees by combining the above formulas: 

So = 7 (2 -W%), (3.16) 

where 



e a = T/ (z,-Waz + ^^( ] ^^-Sj]z 7 ) r (3.17) 



Tf = -fEjlft, W ft = lL a .W ! ', W a = 6 ttfi WP, (3.18) 
W 2 = W a W a = W'Wi = l-\. (3.19) 

The interpretation of these expressions is clear. Each spacetime point is the intersection of 



9 The Zermelo problem is formulated as follows: find the minimal-time navigation road on a geometry d(. 2 = 
hjjdx'dxl, in the presence of a moving fluid creating a drift current (wind or tide) W = W3, with a ship of 
fixed propelling velocity (i.e. fixed with respect to the frame comoving with the fluid or, put differently, of given 
power). The answer is reached by searching for null geodesies of <3.13l l. 
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two lines, belonging each to the two congruences under consideration. At this point W a are 
the spatial velocity components of the inertial observer in the spatial frame of the accelerated 
observer and V7 2 = 1 — VV 2 the corresponding Lorentz factor. 

It is worth making several further comments. The synchronous hypersurface for the 
inertial observer is by definition dual to the time vector. Since df(9,) = (equivalent to 
z°(z a ) = 0), this hypersurface is spanned by {9;}, and hence is not orthogonal to the in- 
ertial congruence (d t • 9 ; - = b\). The orthogonal lines to the Papapetrou-Randers (inertial) 
observer's synchronous hypersurface are nothing but the accelerated congruence tangent to 
the vector field zo (defining the corresponding Zermelo frame) because zq ■ 3, = 0, whereas 
the hyperplanes orthogonal to the the inertial congruence (tangent to eo = 9f) are spanned 
by {bid t + 9;}. Therefore, since (dt — b)(dj + bjd t ) = (equivalent to e°(e a ) = 0), the time 
t of Zermelo observers, i.e. the dual of the hypersurface everywhere tangent to the latter 
hyperplanes, would satisfy dr = dt — b. Such a time cannot be defined since db = 2co ^ 0. 
Put differently, no hypersurface exists tangent to the hyperplanes spanned by {bjd t + 9,} - 
Frobenius theorem. This is a well known manifestation of vorticity. 

The last statement again distinguishes the Papapetrou-Randers and Zermelo observers, 
which are otherwise dual to each other. As for the perception of the rotation, Papapetrou- 
Randers observers feel it through embarked gyroscopes (see d3.10|l ), whereas their inertial 
motion as witnessed by Zermelo observers satisfies 

Vz u = oo$ a 6 a %. (3.20) 

Here u = §q is the velocity of the inertial observers, while cof h are the vorticity components 
as observed in the Zermelo frame: ccfin = L^LpOOij and co^ = VJ K co^. Hence, for the ac- 
celerated observers, the inertial ones are subject to a Coriolis force: Zermelo observers are 
rotating themselves. The velocity vector u = So of the inertial observers undergoes a preces- 
sion around the worldline of a Zermelo observer tangent to zq. The latter being accelerated, 
the variation of u is actually better captured as a Fermi derivative along Zq: 

D, u = (a; z a - z a ( 7 )) S*% + W a z fl (7)2:0, (3.21) 

where 2^(7) = The extra terms result from the rotation of the Zermelo frame and 

contribute to the observed precession of the velocity vector u. 

One can try to tune rotating frames so as to make the perceived angular momentum 
of a given congruence disappear, i.e. make its Fermi derivative vanish with respect to the 
rotating frame. This leads to the so called zero angular momentum frames (ZAMO |25l ). 
In general, Zermelo frames are not ZAMO frames for the Papapetrou-Randers congruences, 
as the Fermi derivative d3.21| | is generically non-zero. It can however be zero under the 
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necessary and sufficient condition 

W'coji = 7 d i7> (3.22) 

since this implies that the combination cOq x — z a ( r y) as well as the coefficient of Zq vanish. 
Equation (|3.22| | carries intrinsic information about the background and can indeed be recast 
as 

C h e Q = 0. (3.23) 

When fulfilled, the Zermelo observers coincide with the locally non-rotating (or ZAMO) 
frames |[25| . Remarkably this occurs for a particular case that will be discussed in our sub- 
sequent developments. 

3.3 Further properties and analogue gravity interpretation 

In the above analysis and particularly in the change of frame from Papapetrou-Randers to 
Zermelo, it has been implicitly assumed that W 2 < 1. The velocity of Papapetrou-Randers 
observers with respect to the Zermelo frame is however dictated by the geometry itself since 
W 2 = bib 1 , and nothing a priori guarantees that b\V < 1 everywhere. There are regions 
in x-space where indeed bib' > 1 bounded by a hypersurface where bib' = 1. The latter 
was called velocity-of-light hypersurface in H2T1 since this is the edge where the Papapetrou- 
Randers observer reaches the speed of light with respect to the Zermelo frame. 

The problem raised here is a manifestation of the global hyperbolicity breakdown. In- 
deed, we have seen that in geometries of the Papapetrou-Randers form (I3.1L constant-f 
surfaces are not everywhere spacelike. The extension of the physical domain accessible to 
the inertial observers moving along u = 9 t is limited to spacelike disks in which b 2 < 1 
holds, bounded by the velocity-of-light surface, where these observers become luminal. 

The breaking of hyperbolicity is usually accompanied with the appearance of closed 
timelike curves (CTCs). These are ordinary spacelike circles, lying in constant-f surfaces, 
which become timelike when these surfaces cease being spacelike, i.e. when b 2 > 1. The 
CTCs at hand differ in nature from those due to compact time (as in the SL(2, R) group 
manifold), and cannot be removed by unwrapping time. They require an excision proce- 
dure for consistently removing, if possible, the b 2 > 1 domain, in order to keep a causally 
safe spacetime. This is comparable to what happens in the case of the three-dimensional 
Banados-Teitelboim-Zanelli black hole [26] - although in the latter case the trouble is not 
due to hyperbolicity issues. We will come back to the CTCs when studying the anti-de Sitter 
Taub-NUT geometry. 

Although the issue of hyperbolicity is intrinsic to our stationary geometries, moving from 
the form d37Q to the form (13131 may provide alternative or complementary views. In the 
Zermelo form (I3.13D the trouble is basically encapsulated in the conformal factor. However, 
some problems such as the original Zermelo navigation problem referred to in footnote [9] 
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are sensitive to the general conformal class^ of (I3.13I I, and the conformal factor j 2 can be 
dropped or replaced. Doing so can leave us with a geometry potentially sensible every- 
where. This instance appears precisely in analogue gravity systems. 

Metrics of the form (I3.13D are in fact known as acoustic or optical (see the original works 
[2711281 or e.g. ITTT1 for an up-to-date review). They are used for describing the propagation 
of sound /light disturbances in relativistic or non-relativistic fluids moving with velocity 
VP in spatial geometry hu, and subject to appropriate thermodynamic/hydrodynamic as- 
sumptions. In this approach, the full metric (I3.13D is an analogue metric and is not the actual 
metric of physical spacetime. Under this perspective, peculiarities such as CTCs, potentially 
present in the analogue geometry, have no real, physical existence. They are manifesta- 
tions of other underlying physical properties such as supersonic/ superluminal regimes in 
the flowing medium. 

In order to be concrete, we would like to quote two examples. The first is the original one, 
where a fluid is flowing on a geometry dl 2 = hjjdx'dxi with velocity field W = VPdj. Assum- 
ing the fluid is non-relativistic, isolated, with mass density q and pressure p, in barotropic 
evolution (i.e. such that the enthalpy variations satisfy dh = d P/o) and with sound veloc- 
ity c s = V VWdp, one finds that irrotational acoustic perturbations propagate along null 
geodesies of the metric 

ds 2 = ^ (-c 2 s dt 2 + (dx l - Wdt) (dx> - Wdt) ) , (3.24) 



c, 



and satisfy the corresponding scalar field equation. The metric (13.241 1 is of the form d3. 13b ■ 
It is however analogue and not the actual spacetime geometry, which is Galilean. This ana- 
logue metric is called acoustic in the case at hand. Similar results are available for light 
propagation, leading to optical geometries (see e.g. |121|29l ). 

Similarly we can quote the case of a relativistic conformal fluid, at rest in the Papapetrou- 
Randers frame of a Papapetrou-Randers geometry d3.1| | - with B = 1. In this case, fluid lines 
are tangent to u = d t , the kinematics is shearless and expansionless with vorticity (I3.9L and 
thus it is non-dissipative. It is also geodesic (see d3.6l l) and Eqs. (12.251 1 imply that e, p and 
T are constant everywhere. The propagation of irrotational perturbations in this set up is 
captured by the following acoustic metric^: 



TD-2 

ds 2 = ^ d = {-(di - hdx 1 ) 2 + (D — l)aij(x)dx l dx^ , (3.25) 

of the Papapetrou-Randers form, which can be recast in the Zermelo form (I3.13I) , following 
63D and rt3l2l 



10 In the Zermelo navigation problem we look for null geodesies. In that framework, going to regions where 
7 2 < means having a drift current faster than what the ship can overcome. 
11 This result is easy to establish, following e.g. similar reasoning as in 1 30 [ . 
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4 Examples in 2 + 1 dimensions 



Examples of Papapetrou-Randers geometries are numerous, possessing diverse properties 
regarding their isometries, their curvature, the regularity of their Randers-Zermelo transfor- 
mation, etc. Their expressions can be simple in Papapetrou-Randers form and complicated 
in the Zermelo representation or vice versa. All this depends in particular on the dimen- 
sion^ We will here focus on a few three-dimensional examples that turn out to emerge 
as holographic duals of exact four-dimensional bulk spacetimes. Although the nature of 
the boundary three-dimensional spacetime regarding Einstein's equations plays little role in 
holography, some underlying intrinsic properties appear to be generic for the backgrounds 
under consideration, and would deserve further investigation. Furthermore, all examples 
below are homogeneous spaces^E even though neither was homogeneity an a priori crite- 
rion, nor did our list provide an exhaustive classification of stationary backgrounds with a 
spatially homogeneous timelike Killing vector field. 

4.1 Warped three-spheres: Bianchi IX 

Warped three-spheres, earlier more appropriately called biaxially squashed three-spheres, 
are deformations of the standard homogeneous and isotropic (round) S 3 . This deformation 
breaks the original SU(2) x SU(2) isometry group down to SU(2) or SU(2) x 11(1). These 
spaces can be endowed with the metric 

d S 2 = f;(7^) 2 ' (4-1) 

1=1 

where 7; are constants and a 1 are the left-invariant Maurer-Cartan forms of SU(2). In terms 
of Euler angles 0<6><7r,0<</>< 271,0 < ip < 4n, these one-forms read: 

a 1 = sin 6 sin xp dcp + cos ip d6 
< a 2 — sin 6 cos ip dcp — sin ip d0 (4.2) 
a 3 = cos 8 dip + dtp. 

For reasons that will become clear in the following, we will consider situations where 
71 = 72- The spaces obtained in this way are called Berger spheres. They are axisymmetric 
i.e. have SU(2) x U(l) isometry group. Since we are interested in spaces with Lorentzian 

12 See [21 J for a detailed account of properties and examples in 3 + 1 dimensions. 

"Following 1311321 , homo geneous three-manifolds include all 9 Bianchi groups plus 3 coset spaces, which are 
H 3 , H 2 x S 1 , S 2 x S 1 (S" and H n are spheres and hyperbolic spaces respectively). 
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signature, we must set negative 7 2 and the metric finally reads: 
ds 2 = L 2 Ua 1 ) 1 + {o 2 ) 2 ] -An 1 (a 3 ) 2 = —An 1 (dtp + cos 9d(p) 2 + L 2 (d9 2 + sin 2 9d<p 2 ) , (4.3) 



where 2L is the radius of the original undeformed S 3 and Ink the deformation parameter 
(k = Vl). 

The time coordinate in (I4.3D is In the original Euclidean sphere this was an angle, 
but there is no reason to keep it compact in the Lorentzian version at hancJ^l Introducing a 
non-compact time t = —2n(xp + cp), the metric (I4.3I I assumes the form: 

ds 2 = - (dt + An sin 2 8/2 d<pf + L 2 (d6 2 + sin 2 0d<£ 2 ) . (4.4) 

This metric is of Papapetrou-Randers type ( I3.lt . The base dx 2 = a„(x)dx'd;^ is a two-sphere 
of radius L, while b = — 4n sin 2 6 / 2 d</> is a Dirac-monopole-like potential. The latter creates 
a constant - in the Papapetrou-Randers orthonormal coframe - vorticity 

co = -nsm9d9 Ad(p = -nk 1 e l Ae 1 (4.5) 

for the geodesic congruence tangent to u = = d t (up to a delta-function contribution of 
the "Misner point" at the southern pole - this name will be justified later). 

As already quoted, the space at hand is homogeneous and belongs to the family of spaces 
invariant under a four-parameter group of motions 1133 U34L here generated by the following 
Killing vectors: 

£1 = — sin</>cot030 + cos(pd0 — 2n^|(l — cos0) d t 
£2 = cos (p cot 9 d,p + sin <p d e + 2n^| (1 - cos 9) d t 
£ 3 = d<p- 2n d t 
e 3 = -2nd t . 

The three former generate the left SIT (2), whereas the latter generates an extra R factor 
(instead of U(l) since t is non-compact). 

The background jA.A) is not globally hyperbolic. Even though it is homogeneous, constant- 
t surfaces are not and 7 = V yJi-in 2 k 2 i<m 2 <>/2 diverges when 9 reaches 9* = 2 arctan L /in. Hy- 
perbolicity holds in the disk < 9 < 9*, whereas it breaks down in the complementary disk 
(0* < 9 < 7r) centered at the Misner point, where 9^ becomes timelike. As a consequence, 
the circles tangent to become CTCs for 9* < 9 < n. Homogeneity implies furthermore 
that CTCs are present everywhere, passing through any arbitrary point of spacetime. In par- 

14 This point deserves some comments, which we postpone to the discussion on the bulk geometries leading 
to (43) as boundary. 
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ticular, for < 8 < 6* the CTCs are sections of cylinders normal to the constant-f surfaces. 
The time coordinate t evolves periodically along these elliptically shaped CTCs. 

As we will see in the forthcoming sections, the situation described here is quite generic 
for three-dimensional homogeneous spacetimes. These include the case of Som-Raychaud- 
huri (Bianchi II) and the celebrated Godel space (Bianchi VIII). They are illustrative examples 
of how homogeneity combined with rotation often leads to the breakdown of hyperbolicity 
and the emergence of CTCs. Godel space in particular was the first to be recognized as 
plagued by CTCs. The CTCs present in these spaces, however, are not geodesies l34ti36| . 
Their presence is therefore harmless for classical causality. This is why Godel-like solutions 
like the case under consideration have never been truly discarded, leaving open the possi- 
bility of quantum mechanical validity^ 

Let us incidentally mention that the above three-dimensional geometry (|4.4|) - when 
lifted to four dimensions by taking the direct product with an extra flat direction - has been 
shown to satisfy Einstein's equations with cosmological constant and energy-momentum 
tensor produced by some specific charged fluid |34| . Alternatively, it also satisfies the equa- 
tions of topologically massive gravity 11421 , subject to a Kerr-Schild |46| deformatior0 

R m 1 fc 2 (l+4n 2 fc 2 ) 
R ¥ v ~ + A^g }lv = ^C ¥V + L UjlUv> (4.7) 

where u = —e^ = — dt — 4n sin 2 Q/idcj), = k2 /4 is the cosmological constant for an unde- 
formed S 3 and C fll , are the components of the Cotton-York tensor defined as 



C f = ^-VJRV^]. (4.8) 



1 

4* 



For the background (|4.4|) this tensor is 

C w dxW = nk A (1 + An 2 k 2 ) [2u 2 + L 2 (d6> 2 + sin 2 6d(p 2 )] . (4.9) 
For later convenience we also quote: 

R w dxW = 2k 4 n 2 H 2 + (l + 2n 2 k 2 ) (d9 2 + sin 2 Odcp 2 ) (4.10) 

and 

R = 2k 2 (1 + n 2 k 2 ) . (4.11) 



15 Attempts, among others in string theory within holography, were proposed a few years ago (see e.g. I36H4T1 
and references therein). 

16 It should be mentioned that the metric <4.4t solves also topological massive gravity equations without Kerr- 
Schild deformation (43lf45l , provided one tunes appropriately the relationship among L 2 and k 2 . This holds for 
generic warped homogeneous spaces like AdS3, studied in Sec. 14.21 
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4.2 Warped AdS 3 : Bianchi VIII 

Following the paradigm of squashed three-spheres, studied in Sec. 14.11 we analyze here the 
deformations of AdS3. The latter is the (universal covering of the) SL(2, R) group manifold 
and has left and right SL(2,R) isometries. Homogeneous deformations can break partially 
or completely one SL(2,R) factor. 

There are many realizations of the Maurer-Cartan forms of SL(2, R). The one we choose 
here is convenient for the specific deformation we will consider in the following: 

p° = — dr + cosher dep 
< p 1 = — sin t da — sinh a cos T d(p (4-12) 
p 2 = cos t da + sinh a sin T dep, 

where < (p < 2n, < c < +oo, and T G [0,271] or better R if we consider the universal 
covering of the space. The metrics under consideration are of the form 

ds 2 = - ( 7o p ) 2 + (71P 1 ) 2 + (72P 2 ) 2 • (4.13) 

The 7s being constant, this geometry is homogeneous. When V/, 7; = L, we recover radius- 
2L AdS 3 . 

We will restrict ourselves here to the elliptically-squashed AdS3, obtained with 70 = 2p S 
R and 71 = 72 = L. These geometries have a 4-parameter isometry group SL(2, R) x R. The 
Abelian factor is the remaining one-parameter subgroup of the broken isotropy symmetry - 
(J(l) before taking the universal covering. 

Elliptic deformations were studied as such in [47] and the continuous line obtained there 
coincides - when uplifted to four dimensions - with the family of spacetime-homogeneous 
Godel-type metrics discussed in |34| . After a coordinate transformation trading x for t = 
2p(r - <p) € R, the metric (14.131 1 reads: 

ds 2 = - (dt - 4p sinh 2 <r/% d(p\ 2 + L 2 (da 2 + sinh 2 a d<£ 2 ) . (4.14) 

This is a timelike fibration over a hyperbolic plane H2. 

Alternative inequivalent deformations exist in SL(2,R), such as the hyperbolic or the 
parabolic ones which lead to spacelike fibrations over AdS2 or plane-wave superpositions 
with AdS3. They have been studied extensively in the framework of string theory [|38ti40| 
and more recently in holography [|45| . However, these are not of the Papapetrou-Randers 
form d3.lt , as opposed to the elliptic deformation for which b = 4p sinh 2 a /i d(p while the 
base dx 2 = ajj(x)dx'dxi is radius-L Lobatchevsky plane. The "non-compact Dirac potential" 
b creates for the geodesic congruence tangent to u = Sq = d t a homogeneous vorticity which 
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- in the Papapetrou-Randers orthonormal coframe - reads: 



oj = p sinh a da A d<p = pk 2 e 1 A e 2 . 



(4.15) 



In the case at hand, the base is non-compact and the Misner point is rejected to infinity. 

The hyperbolicity properties of the background (|4.14|l are richer than for (|4.4[) . Indeed, for 
the space at hand, the Lorentz factor relating Zermelo and Papapetrou-Randers observers is 
7 = V v / i-4p 2 fc 2 tanh 2 <V2. This factor remains finite for finite a provided p < L /2- the limiting 
case p = L /i corresponding to the undeformed AdS3. Under this condition, the space is 
globally hyperbolic and constants surfaces are spacelike everywhere. When p > L /i this 
property breaks down for a > a* = 2arctanh V 2p. Constant-f surfaces are spacelike on 
a disk only, centered at a = and bounded by a velocity-of-light surface located at a = 
(T* . The breakdown of hyperbolicity is accompanied with the appearance of CTCs, present 
everywhere as a consequence of homogeneity. This happens in particular for p = L / V2 
corresponding to Godel's solution. 

As mentioned earlier, the absence of hyperbolicity in the backgrounds considered in our 
general framework is closely related to the combination of vorticity and homogeneity. This 
is even better illustrated in the AdS3, where the non-compact nature of the base makes it 
possible to evade the breakdown, provided the vorticity is small enough with respect to 
the scale L set by the curvature. The isometry group is in the present case SL(2, R) x R 
generated by the following Killing vectors: 



£o = d<p ~ 2p 9 t 

l\ = - cos <p coth cr d<p - sin <p d a - 2p^^ (1 - cosh a) d t 
fr = - sin (p coth 0-3$ + cos (pd a - 2p4^ (1 - cosh a) d t 
qo = -2pd t . 



(4.16) 



For completeness we also quote the Levi-Civita curvature properties of the metric (|4.14[) : 



C pv dxW = pk 4 (1 - 4p 2 fc 2 ) 2u 2 + L 2 (do 2 + sinh 2 a d(p 2 
R ¥V dx''dx" = 2k 4 p 2 u 2 - (1 - 2p 2 it 2 ) (da 2 + sinh 2 adcp 2 ^} , 
-2k 2 (1 - p 2 k 2 ) , 



R 



(4.17) 

(4.18) 
(4.19) 



where u = — dt + 4p sinh 11 1 2 dep. 



4.3 Rotating Einstein universes 

In the limit n — > 0, the geometry (|4.4|l becomes a homogeneous metric on R x S 2 . This is 
the Einstein static universe, which is a trivial case of static Papapetrou-Randers geometry. 
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The isometry group remains unaltered and generated by the Killing vectors (I4.6I I. The metric 
reads: 

ds 2 = -dt 2 + L 2 (d9 2 + sin 2 9 d<p 2 ) . (4.20) 
Trading (9,(p) for (6',(p') defined as 



</> = <// + Cl^t, D e A , = E, 2 = 1- L 2 nl, 
D e = \- L 2 Q 2 sin 2 Q, A e , = 1 - L 2 f& cos 2 & , 



(4.21) 



we obtain 



ds 2 = -dt 2 + ^ ^df? /2 + sin 2 0' [d<// + Qeodf] ^ . (4.22) 



This is the Einstein universe uniformly rotating at angular velocity Ooo in spheroidal coordi- 
nates. 

The metric (|4.22|l is in a Zermelo form d3.13t . It can be brought to Papapetrou-Randers 
form using the general transformations (13.111 1 and (I3.12D : 



ds 



2 _ 

As/ 



dt 



L 2 Q 



sin d</> 



L 2 d9' 2 L 2 A#/ . 2 , (2 



As/ 



sin z 0'df " 



(4.23) 



In order to reach a canonical Papapetrou-Randers form d3.1| | with B = 1, a conformal trans- 
formation is required with conformal factor <3> = A s'/ a. The resulting geometry 



ds 2 



dt 



L 2 Q r 



sin 2 6'd(p' 



+ 



L 2 d0 



12 



L 2 A a , 



sin 2 e'd(p' 2 , 



(4.24) 



appears as the boundary of Kerr-AdS four-dimensional bulk geometry [48 J. We would like 
to stress that although (I4.23I I is still invariant under SU(2) x R generated by (I4.6I I (at n = 0), 
this is no longer true for (|4.24|l , where the isometry group is reduced to U(l) x R generated 
by d<p and d t , while £i and £2 in (I4.6D become conformal Killing vectors. 

Let us finally mention that we can in a similar fashion consider the limit p — > in the el- 
liptically deformed AdS3 geometries, given in (|4.14[) . This leads to a Einstein-static-universe- 
like geometry of the type R x H2, which can be brought into the form 



ds 2 



Z 



dt 



L 2 fL 



sinh 2 (fd(p' 



iJxkr' 

©n 



a 



^sinhVd^' 2 



(4.25) 



after a coordinate transformation (o~,(p) h4 (o~',(p') with 

<p = <p' + fW, Ho-Qpi = Z, Z = 1 + L 2 Q 2 C/ 
H a = 1 — L 2 2 sinh 2 cr, 0^ = 1 + L 2 Q 2 <J cosh 2 (/. 



(4.26) 
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One can similarly perform a conformal transformation leading to 



ds 2 



df 



L 2 O c 



sinh a' dtp' 



L 2 dcr 



©n 



2 +^ sinh Vdtf/ 2 . 



(4.27) 



The geometry described by (|4.27|l is globally hyperbolic without any restriction. For the 
metric (|4.24[) , global hyperbolicity holds under the condition O,*, < k. Again this is a bound 
on the magnitude of the vorticity 



L 2 Q co sin20 / , , 
co = — d6 A dip 



(4.28) 



carried by the geodesic congruence u = Sq = d t in this background - and is relaxed for the 
non-compact base (squashed H2 in (|4.27| l versus squashed S 2 in (|4.24|l ), where 



L 2 Q 0O srnh2cr / , , 

co = — da A dtp . 

2Z r 



(4.29) 



4.4 Zooming at the poles: Bianchi II and VIIq 

We can further analyze the geometries met so far by zooming around their poles. This ex- 
hibits new backgrounds of Papapetrou-Randers type, interesting for their own right. 



Som-Raychaudhuri and Heisenberg algebra. The geometry around 6 in (|4.4|l or a 
in (|4.14|l is 

ds 2 w - (dt - L 2 n 0O ^ 2 d</>) 2 + L 2 (d^ 2 + x 2 dcp 2 ) (4.30) 

with x — $ or a, and Qoo = —nk 2 or pk 2 . This very same geometry appears also around 
9 « or 7T in ()4.24|l and around (7 ~ in (|4.27|l - upon appropriate definitions of x involving 
E or Z. 

Metric d4.3QI > is the Som-Raychaudhuri space, found in [49 1 by solving Einstein equations 
with rotating, charged dust with zero Lorentz force. It belongs to the general family of three- 
dimensional homogeneous spaces possessing 4 isometries studied in [33,34], which include 
the various metrics that we have discussed so far here. In the case of Som-Raychaudhuri 
(Eq. (14.30 l l) the isometries are generated by the following Killing vectors: 

K x = k^£ dcf, — k cos tpd x — LOoo^ sin tp dt 

K lJ = k^-d ( p + ksincpd x -LCl co xcos(pd t 
< (4.31) 

Ko = 20.00 dt 
K = d f . 

The vectors K x , K y and Kq form a Heisenberg algebra, and indeed the Som-Raychaudhuri 
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metric can be built as the group manifold of the Heisenberg groupfBianchi II) at an extended- 
symmetry (isotropy) point with an extra symmetry generator 9^. Actually this corre- 
sponds to a contraction of SU(2) x R into a semi-direct product of the Heisenberg group 
with an extra LT(1) generated by K x = —k^\,K y = k^ 2 ,K = k 2 e 3/ K = £ 3 - e 3 (see <|46l l). It 
can similarly emerge from the SL(2,R) x R algebra (I4.16I I. 

Similarly to Godel space, Som-Raychaudhuri space contains non-geodesic closed time- 
like curves. These are circles of radius x larger than 1 /n co |50l . 



Flat vortex and E(2) geometry. The southern pole of (|4.4|) is not captured in the above 
scheme. Indeed this is a fixed point of the transformation generated by the Killing vector 
£3 + £3 (see d4.6l l), which we referred to as the Misner point. Around this point 



ds 2 



[dt + n(4- X 2 )d<p) + L 2 (d* 2 + x 2 d<p 2 ) , 



(4.32) 



where x = tc — Q. The latter is known as a flat vortex geometr y h omogeneous and invariant 
under an E(2) x R algebra (E(2) is Bianchi VIIq) generated bycl 



U = = (fy - 4n d t ) - k cos <p d x 
Uj = -fcfc = ^ (3* - 4n d t ) + k sin <p d A 
L = ^3 = 9<f - 2n d t 



(4.33) 



L — £3 = — 2n d t . 
It also appears as a contraction of the SU(2) x R. 



5 Holography in the Fefferman-Graham expansion: from 3 + 1 to 
2 + 1 

In this section we study the properties of holographic fluids, as they arise in the Fefferman- 
Graham expansion along the holographic radial coordinate. The set of data reached in this 
manner contains the boundary frame i.e. the geometrical background hosting the fluid, as 
well as the energy-momentum tensor, which describes the fluid dynamics. This general 
method is applied to selected four-dimensional solutions of Einstein's equations, whose 
boundaries coincide with Papapetrou-Randers geometries studied in Sec. 0] 

17 [K x ,Ky\ = Ko, [K Xt Ko] = [Ky, Kq] = 0, [K,K X ] = K y , [K, Ky] = -K x and [K,K ] = 0. 
18 [L x , Ly] = 0, [Lq, L x ] = Ly, [Lq, Ly] = — L x and everything commutes with L. 
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5.1 Split formalism and Fefferman-Graham in a nutshell 

We find illuminating to discuss holographic fluid dynamics in D = 2 + 1 dimensions starting 
from the 3 + 1-split formalism introduced in l5T1 - [53"ll . We begin with the Einstein-Hilbert 
action in the Palatini first-order formulation 

r e ABCD (r ab -^E a AE b )aE c AE d 



32tzG n 

where Gn is Newton's constant. We also assume negative cosmological constant expressed 
as A = —3/ l 2 = —3k 2 . We denote the orthonormal coframe E A , A = r,a and use for the bulk 

metric the signature H 1 — h The first direction r is the holographic one and x = (£, x 1 , x 2 ) = 

(t,x). 

Bulk solutions are taken in the Fefferman-Graham form 

ds 2 = ^dr 2 + ^ nab E a {r,x)E b {r,x) . (5.2) 



For torsionless connections there is always a suitable gauge choice such that the metrics (I5.2D 
are fully determined by two coefficients e a and f" in the expansion of the coframe one-forms 
E"(r, x) along the holographic coordinate r 6 K+ 



E a (r,x) 



t(x) + ^F a (x) + - 



L? r 

r 3 



/ fl (x) + --- . (5.3) 



The asymptotic boundary is at r — > oo and it is endowed with the geometry 

ds 2 

d ^=Jfej^a- (5 - 4) 

The ellipses in (I5.3I I denote terms that are multiplied by higher negative powers of r. Their 
coefficients are determined by e a and f a , and have specific geometrical interpretations^, 
though this is not relevant for our discussion. 

The 3 + 1-split formalism makes clear that e a (x) and f a (x), being themselves vector- 
valued one-forms in the boundary are the proper canonical variables playing the role of 
boundary "coordinate" and "momentum" for the (hyperbolic) Hamiltonian evolution along 
r. For the stationary backgrounds under consideration, describing thermally equilibrated 
non-dissipating boundary fluid configurations, e a and f" are f-independent. 

The boundary "coordinate" is given by the set of one-forms e a . For this coframe we must 
determine the "momentum" of the boundary data. For example, when the boundary data 
carry zero mass, we expect this to be zero. In this case f a (x) = and the unique exact 
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For example, the coefficient F a is related to the boundary Schouten tensor. 
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solution of the Einstein's equations is pure ACIS4. 

More generally, the vector-valued one-form f satisfies 

f a Ae a = 0, e abc f a Ae b Ae c = 0, e abc Df h A e c = , (5.5) 

where the action of the generalized exterior derivative D on a vector-valued one-form V" is 
defined as 

DV a = dV a + e\ c B h A e c , (5.6) 

and the "magnetic field" B a is the Levi-Civita spin connection associated with e a |52| . One 
can easily see that conditions (15.51 1 imply respectively symmetry absence of trace and co- 
variant conservation of the tensor T = T a b e a <g> e b , defined as 

/« = lT(e«) = ir/, K = TT^T- ^ 
k k 8tcGnL 

Hence we can interpret the latter as the covariantly conserved energy-momentum tensor of 
a conformal field theory. Here we are interested in particular stationary bulk solutions for 
which we expect the energy-momentum tensor be reduced to the perfect relativistic form, 
Eqs. d2T9l and flZ20} : 

T a b = (e + p)u a u b + pS a b . (5.8) 

Although the two necessary ingredients for the description of a relativistic perfect fluid, 
namely the boundary frame and the velocity one-form, are nicely packaged in the leading 
and subleading independent boundary data, until now we did not assume any specific rela- 
tionship between them. Nevertheless it is clear that such a relationship would be imposed 
by any exact solution of the bulk gravitational equations, given the interior boundary condi- 
tions. We will soon observe that the Fefferman-Graham expansion of the exact solutions of 
Sec. I5.2| yield the same form for the boundary energy-momentum tensor, namely 

The boundary frame one-forms e a are themselves, of course, different in the three solutions. 
Comparing (15.71 1, (15.81 1 and (I5.9L we find 

M 

£ = 2p = 2K— , (5.10) 

constant as already advertised. The solutions under consideration describe thus the same 
conformal fluid in different kinematic states. More importantly, (I5.9I I fixes the direction of the 
velocity field with respect to the boundary frame to be 

u = eg . (5.H) 
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As explained in Sec. 13.11 in the Papapetrou-Randers geometry (Eq. (|3.1| |), this congruence is 
tangent to a constant-norm Killing field and has thus zero shear, expansion and acceleration 
(consistent, according to (I2.25D , with the constant pressure found in (I5.10I I). It also shows 
that the observer's frame e a is comoving. Therefore, in the Fefferman-Graham expansion the 
kinematic properties of holographic fluids are determined by the geometric properties of the 
boundary comoving frame. 

5.2 Some exact geometries and their fluid interpretation 

We present here three examples of holographic fluids with vorticity that reside at the bound- 
ary of two exact solutions of the bulk vacuum Einstein equations; the Kerr-AdS4, the Taub- 
NUT-AdS 4 and the hyperbolic NUT-AdS 4 black hole solutions. 

The four-dimensional Kerr solution of Einstein's equation with A = —3k 2 reads: 



ds 2 



dr 2 



V{r,6) 



V{r,6) dt--sin 2 6d(p 



12 p2 



+ ^-de 2 



sin 2 9 At 



adt 



r 2 + a 2 



dtp 



,(5.12) 



where 



and 



A,, = (r 2 + a 2 )(1 + k 2 r 2 ) - 2Mr 



p 2 = r 2 + a 2 cos 2 1 



A e = 1 - k 2 a 2 cos 2 9 
S = 1 — k 2 a 2 . 



(5.13) 



(5.14) 



The geometry has inner (r_) and outer (r+) horizons, where A,- vanishes, as well as an ergo- 
sphere at g# = 0. The solution at hand describes the field generated by a mass M rotating 
with an angular velocity 

a(l +k 2 r\) 



r\ +a 2 



(5.15) 



as measured by a static observer at infinity Il54ll55ll . Note that asymptotic observer associ- 
ated with a natural frame of the coordinate system at hand is not static, but has an angular 
velocity 

Qoo = ak 2 . (5.16) 



The boundary metric of Eq. (I5.12D is (|4.24[) (without primes in the angular coordinates). 

The Taub-NUT-AdS4 geometry is a foliation over squashed three-spheres solving Ein- 
stein's equations with negative cosmological constant (the as are given in (|4.2| |): 



ds 2 



dr 2 

W) 
dr 2 
V(r) 



(^) 2 + (^ 2 



An 2 V(r) (a 3 )' 



+ (r 2 + n 2 ) (d9 2 + sin 2 9 d(p 2 ) - V(r) 



dt + 4n sin 2 - d(p 



(5.17) 
(5.18) 



23 



where 



(5.19) 



and 



A,. = r 2 - n 2 - 2Mr + k 2 (r 4 + 6nV - 3n 4 ) 
p 2 = r 2 + w 2 . 



(5.20) 



Besides the mass M and the cosmological constant A = —3k 2 , this solution depends on an 
extra parameter n: the nut charge. 

The solution at hand has generically two horizons (V(r±) = 0) and is well-defined out- 
side the outer horizon r+, where V(r) > 0. The nut is the endpoint of a Misner string H56l , 
departing from r = r + , all the way to r — > oo, on the southern pole at 6 = n. The geometry 
is nowhere singular along the Misner string, which appears as a coordinate artifact much 
like the Dirac string of a magnetic monopole is a gauge artifact. In order for this string to be 
invisible, coordinate transformations displacing the string must be univalued everywhere, 
which is achieved by requiring the periodicity condition t = t — 8nn. 

One can avoid periodic time and keep the Misner string as part of the geometry. This 
semi-infinite spike appears then as a source of angular momentum, integrating to zero |57l 
|58l , and movable at wish using the transformations generated by the above vectors. This 
will be our viewpoint throughout this work. However, despite the non-compact time, the 
Taub-NUT-AdS geometry is plagued with closed timelike curves, which disappear only 
in the vacuum limit k — > 1591 . Even though this is usually an unwanted situation, it is 
not sufficient for rejecting the geometry, which from the holographic perspective has many 
interesting and novel features. The boundary geometry of (I5.18L reached by following (I5.4L 
is (|4.4|l and this is the background where the fluid is evolving. What we called Misner point 
for the latter in Sec. l4.ll is the endpoint of the bulk Misner string. 

On a non-compact horizon, the nut charge can be pushed to infinity. This happens in 
hyperbolic NUT black holes, obtained as foliations over three-dimensional anti-de Sitter 
spacea§§. Using the SL(2,R) Maurer-Cartan forms (|4.12[) , we obtain the following solution 
of Einstein's equations with cosmological constant A = —3k 2 : 



20 In their Euclidean section these geometries have no nut, but only a bolt (60"]. We shall nevertheless conform 
to standard use and call them - with a slight abuse of language - hyperbolic NUT black holes, to stress the presence 
of a non-trivial S 1 fibration over Hz- Physically, they represent rotating black hyperboloid membranes |61 1. 




V(r) 
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with V(r) given in (I5.19I I and 



= -r 2 + p 2 - 2Mr + k 2 (r 4 + 6p 2 r 2 - 3p 4 ) 
= r + p . 

Here M is the mass parameter and p characterizes the non-trivial S 1 fibration over the H2 
base. In this case no Misner string is however present, and the space is globally hyperbolic 
provided p < L /i. The boundary metric of (15.221 1 is 1(4.14)1 , which plays the role of host for the 
holographic fluid. Interestingly, this family of solutions is connected to the Kerr-AdS4 black 
hole, as we will now show. 

The Kerr-AdS4 black hole (I5.12I I has a rotation parameter a restricted to a 2 < L 2 , and is 
singular for a 2 = L 2 . It has however a finite, maximally spinning limit if the a — > L limit is 
taken keeping the horizon size finite and simultaneously zooming into the pole |6l| . More 
explicitly, we trade the angle 9 in (I5.12I I for a new coordinate a according to 



sm 



9 = VM, sinh 2 0/2. (5.24) 



Then, the resulting metric has a regular a — > L limit, 

ds2 = W) + ( dal + sinh2 ^ ~ v{r) 



1 • , 2^ , I 2 

at — — sinh — ad> 
k 2 Y 



(5.25) 



with 



"W-I + A'-ot- (5-26) 

This is again a rotating black hyperboloid membrane. Indeed, performing the rescaling 
r 1— > 2r, t \— > t/ 2, this metric is cast in the form (|5.22|l - (|5.23|l with parameters M = M /8 and 
p = L /2, and the boundary of (I5.25D takes the form (|4.14[) with this value of p. The boundary 
is therefore the undeformed AdS3 spacetime that can equivalently be recovered by perform- 
ing the coordinate transformation d5.24|) directly on the Kerr-AdS boundary metric (|4.24|l , 
followed by the a —¥ L limit. As we saw, this special value corresponds to the limiting 
case that enjoys global hyperbolicity. Accordingly, from the Kerr-AdS point of view, it is 
obtained as the ultraspinning limit for which the boundary Einstein static universe rotates 
effectively at the speed of light |48ll . In this way, we nicely connected fluids on warped AdS3 
backgrounds to fluids living on a rotating Einstein static universe into a single, continuous 
family. 

There are many other exact bulk four-dimensional geometries that one could study under 
the perspective of describing boundary fluids with vorticity. One can for instance consider 
the case of flat horizon reached when the trigonometric sinus in (I5.18I I (or the hyperbolic 
sinus in (I5.22D ) is traded for a linear function and the potential adapted by dropping the first 



25 



two terms (see e.g. l59ll60H ). This bulk solution, also plagued by the global hyperbolicity 
problem, gives rise to a boundary fluid moving on the Som-Raychaudhuri geometry d4.301 > . 
Alternatively, the flat-horizon bulk solution can be obtained as an appropriate pole-zooming 
of the four-dimensional Kerr-AdS, consistent with the observed relationships among the 
boundary geometries. Limits at n, p — > lead to the so-called topological black holes I62ti64| . 
These are interesting in their own right |65] even though the holographic fluid dynamics has 
no intrinsic vorticity - their boundaries are Einstein static universes KxS 2 given in (|4.20| | or 
R x H.2, where no fiber appears that would create vorticity. Hyperbolic Kerr-AdS or other 
exact bulk metrics can be found to reproduce on the boundary (|4.27| l or d4.32| | Il66| . One 
can also find solutions that combine nut charge and ordinary rotation fl67| such as Kerr- 
Taub-NUT-AdS. We will neither pursue any longer the general analysis of this rich web of 
backgrounds exhibiting many interrelations, nor delve into a quantitative presentation, but 
move instead into another interesting approach to holographic fluids, which can be easily 
exemplified with the backgrounds at hand. 

The above bulk geometries describe holographically a conformal fluid at rest without 
shear and expansion in the Papapetrou-Randers frame of a Papapetrou-Randers geometry 
d3.1| |. These fluids have a non-trivial kinematics, though, because of the vorticity of the 
geodesic congruence they fill. The vorticity is different in the various cases: it is given in 
Eqs. (14281 d45t and d415l for Kerr-AdS, Taub-NUT-AdS and hyperbolic NUT-AdS. In 
the first case, the fluid undergoes a cyclonic motion with maximal vorticity at the poles 
and vanishing at the equator. In the other two backgrounds, the vorticity is constant as a 
consequence of the homogeneity. The velocity fields are not homogeneous, though, and 
behave differently in Taub-NUT-AdS and hyperbolic NUT-AdS. 

Even though the boundary spacetime of Taub-NUT-AdS is homogeneous, the constant-f 
surfaces are not. Inertial observers, comoving with the fluid have therefore a different per- 
ception depending on whether they are at < 8 < 8* or in the disk 6* < 8 < n, surrounding 
the Misner string. This gives a physical existence to the b 2 = 1 edge, the meaning of which 
is better expressed in the Zermelo frame. In the latter, the fluid becomes superluminal and 
the Misner string is interpreted as the core of the vortex with homogeneous vorticity. 

The various troublesome features which appear in Godel-like spaces as the ones at hand 
are intimately related with the non-trivial rotational properties combined with the homo- 
geneous character of these manifolds. In other words, for the Taub-NUT-AdS boundary, 
they are due to the existence of a monopole-like Misner vortex^! Although no satisfactory 
physical meaning has ever been given to Godel-like spaces, the causal consistency of the 
latter being still questionable, they seem from our holographic perspective to admit a sensi- 

21 Since the bulk theory is such that the boundary does not have access to a charge current, the Misner vortex 
cannot be associated with a vortex in an ordinary superfluid, but is related to the spinning string of |68|, the 
metric of which, Eq. (432), indeed appears when zooming in on the southern pole. 
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ble interpretation in terms of conformal fluids evolving in homogeneous vortices (14.511 I or 
(OBI) . 

The case of hyperbolic NUT-AdS, Eq. (I5.22D , is yet of a different nature. This bulk ge- 
ometry leads again to homogeneous boundary (|4.14| l. Hence, the fluid has constant vorticity 
(I4.15D . However, in the case at hand, the spatial sections dx 2 = a,,dx ! dx^ are non-compact 
and negatively curved as opposed to the boundary of Taub-NUT-AdS. As a consequence, 
the combination of vorticity and homogeneity does not break global hyperbolicity as long as 
p < L /i. In this regime, the velocity of the fluid is well-defined everywhere, and its Lorentz 
factor with respect to Zermelo observer is increasing with a and bounded as 7 < V ^i-ip 2 k 2 . 
For this observer, the fluid is at rest in the center (i.e. at the north pole) and fast rotating at 
infinity. When p = L /i, it reaches the speed of light when a — > 00, whereas for p > L /i this 
happens at finite a = a*, along the surface-of-light edge. The latter situation is similar to 
what happens in the Taub-NUT-AdS irrespectively of the value of the nut charge n. In the 
hyperbolic case, the major differerence is however that the vortex, together with the Misner 
point are sent to spatial infinity (a — > 00). 

The above discussion holds in the perspective of interpreting the holographic data as 
a genuine stationary fluid. There is however an alternative viewpoint already advertised, 
consisting in the analogue gravity interpretation of the boundary gravitational background. 
From the latter, the physical data are still (hjj, W) i.e. a two-dimensional geometry and a 
velocity field. However, their combination into (13.131 1 is not a physical spacetime. The 
would be light cone, in particular, is narrowed down to the sound or light velocities in the 
medium under consideration - necessarily smaller than the velocity of light in vacuum. Con- 
sequently, the breaking of hyperbolicity or the appearance of CTCs are not issues of concern, 
and the regions where 7 becomes imaginary keep having a satisfactory physical interpre- 
tation as portions of space, where the medium is supersonic/ superluminal with respect to 
the sound/light velocity in the medium and not in the vacuum. Finally, the virtual spacetime 
(13.131 1 governs the mode propagation through the fluid. This way of thinking opens up a 
new chapter that requires adjusting suitably the standard holographic dictionary. The latter 
provides indirect information on the physical system that must be retrieved. 

6 Alternative expansion: from 2 + 1 to 3 + 1 

The holographic fluids we have described so far emerge from exact bulk four-dimensional 
solutions of vacuum Einstein's equations with negative cosmological constant. They appear 
as a set of two pieces of boundary data - the coframe and the energy-momentum current 
- following the Fefferman-Graham expansion at large (appropriately chosen) radial coordi- 
nate. This is a top-down approach as opposed to the alternative bottom-up method initiated 

22 As already stressed, one should add a ^-function contribution to the Taub-NUT-AdS vorticity 44.5t because 
we keep the Misner string physical with non-compact time 1 221231- 
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in |2j. The latter aims at reconstructing perturbatively a bulk solution starting from bound- 
ary data. The perturbative expansion of the bulk geometry obtained in this way is, however, 
orthogonal in spirit to that of the Fefferman-Graham series, since it is a derivative (of the 
velocity field) rather than a large-radius expansion. It captures therefore from the very first 
order the presence of a regular horizon of the black object that generates the dynamics of the 
boundary fluid. 

As a matter of principle, one could generally follow the above procedure and pertur- 
batively reconstruct the bulk solution corresponding to any boundary background of the 
Papapetrou-Randers type in terms of the data containing the full dynamics of the 

fluid. Our viewpoint is different though, and we are here interested in modestly discussing 
the interplay between the perturbative expansion developed in |2l[16j[69l and the exact so- 
lutions we have analyzed in Sec. [5] This is motivated by the observation made in Ifl6ll69l , 
according to which the proposed bulk perturbative reconstruction of the Kerr-AdS bound- 
ary fluid (in several dimensions) does coincide exactly with the original bulk geometry at 
second order - modulo a specific resummation, indicative of the genuinely all-order nature 
of the solution. This is remarkable and leads to the deeper question: why and under which 
conditions does this occur? 

The question raised here is twofold. Given an exact bulk solution, what can make it be 
expressed in the form of a limited expansion in terms of its boundary data obtained via the 
Fefferman-Graham procedure? Given a set of arbitrary boundary data, what could ensure 
the corresponding bulk series be exact at finite order? 

Making progress in this direction would require delving into the physics of dissipative 
phenomena and their holographic expressions. This analysis stands beyond our present 
scope. We can nevertheless make an observation that might ultimately be relevant. In all 
backgrounds under consideration, the bulk can be expressed exactly as a limited deriva- 
tive expansion provided an extra term (with respect to the expansion proposed in |[T6ll69l ) 
involving the Cotton tensor of the boundary geometry is appropriately added, and after per- 
forming a resummation similar to that of Kerr-AdS. This holds for all Papapetrou-Randers 
backgrounds presented in Sec. HI To keep our presentation compact, we will only consider 
those for which we studied the bulk realisation in Sec. |5j namely Kerr-AdS, Taub-NUT-AdS 
and hyperbolic NUT-AdS. 

The starting point for this analysis is the expression in Eddington-Finkelstein coordinates 
of the bulk metrics. For Kerr-AdS, Eq. (I5.12I I this is achieved by performing the following 
coordinate change: 

df i-> dt - 1 ^dr 

L (6-1) 
d(p i-> d(p — ^dr 

with all quantities defined in d520l) . Similarly for Taub-NUT-AdS, Eqs. d5l8l and d5T20T >, one 
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performs 

r 2 + n 2 

dt^dt — dr, (6.2) 

A r 

while the same holds with n H- p for the hyperbolic NUT-AdS given in (|5.22|) and (|5.23|) . All 
three bulk metrics assume then the following generic form: 

ds 2 = -2udr + r 2 k 2 dsl ry + ^ L^dxW + — f 2Mr + ^fe^ J ' (63 > 

where all the quantities refer to the boundary metric ds^ ry . The latter being of the Papapetrou- 
Randers form d3.lL u = — dt + b and to = \db. Furthermore are the components of the 
Cotton tensor (zero for Kerr-AdS, and displayed in Eqs. (I4.9D and (|4.17|l for the other cases). 
Finally 

Z^dxW = -2uV v co v }l dxV - (v/co Xv dxW - u 2 ^, (6.4) 
P 2 = r 2 + ^w 1lv cv^ r (6.5) 

and p 2 , as computed in (I6.5D , coincides with the quantities defined in Eqs. (I5.14D , (I5.20D and 
(15231 . 

The above result (I6.3D deserves a discussion. It appears as a limited derivative expan- 
sion on the velocity field of the boundary geodesic congruence u = d t . The latter carries 
neither expansion, nor shear (see Sees. 12.11 and 13.11 1 but only vorticity given in Eqs. (14.281 1, 
(I4.5I I and (I4.15I I for the three backgrounds under investigation Kerr-AdS (I5.12L Taub-NUT- 
AdS (15.181) and the hyperbolic NUT-AdS d5.22| |. It seems that at most two derivatives of 
the velocity field are involved, but this counting is naive. Indeed, the vorticity being ulti- 
mately an intrinsic property of the geometry^!, R ~ cv 2 , while Ceo ~ a; 4 . Furthermore, 1 /p 2 
is a resummed power series in even powers of the vorticity. As already advertised, this re- 
summation betrays the infinite perturbative expansion underlying the method, that would 
otherwise appear as limited to the fourth order. 

The metric (16.31 1 yields an exact solution of AdS4 gravity for a large class of boundary 
Randers data In addition to the cases described above, one can rewrite in this form 

the full Kerr-Taub-NUT-AdS4 family of metrics, as well as all rotating topological black 
holes found in [|66| : the rotating black cylinder and the rotating hyperbolic black membrane. 
All these metrics belong to the Plebahski-Demiahski type-D class of solutions [|70| . It is an 
interesting problem to find to see if it is possible to extend this collection and find the most 
general set of Randers data (fc„ a,,-) generating an exact solution through (I6.3I I. 

It is remarkable that all known exact AdS4 black hole solutions can be set in the above 

23 Vorticity components are directly related to the connection coefficients, as e.g. T' • = -co'-. 



29 



form (I6.3L much like Kerr-AdS, provided an extra term (one should say an extra resummed 
series) based on the Cotton tensor is added. This term, of fourth order in the derivatives 
of the velocity fielcP^l, was absent in the original expressions of [1611691 , only valid up to 
second order. In five or higher bulk dimensions, terms involving the Weyl tensor appear 
at the second order HUEU/ but obviously do not contribute in the four-dimensional case 
under consideration. Our expression (I6.3I I shows that in four dimensions analogous terms, 
involving the Cotton tensor, appear in the derivative expansion starting from the fourth 
order. 



7 Conclusions 

In this review we presented an extensive discussion of the holographic description of vor- 
ticity. This is the first step in efforts to extend AdS/ CMT to systems such a rotating atomic 
gases of analogue gravity systems. The upshot is that even the simplest setup, namely non- 
dissipating fluids in local equilibrium with non-zero vorticity has an extremely rich geo- 
metric structure whose detailed analysis should lead to new and interesting physical results. 
One such result, presented in 11231 was the calculation of the classical rotational Hall viscosity 
coefficient of neutral 2 + 1 dimensional fluids having uniform vorticity Q 

= ^, (7.1, 

which we were not able to find in recent works on parity broken hydrodynamics in 2 + 1 
dimensions. 

The next steps in our program will certainly reveal interesting physical consequences. 
For example, the study of scalar, vector and ultimately graviton fluctuations around the 
above geometries should lead to the determination of various transport coefficients for ro- 
tating neutral fluids. These developments might also shed light on the thermalization pro- 
cesses that are expected near analogue horizons. Furthermore, we believe that our approach 
offers a well-defined path to study the issue of time-dependence in conjunction with irre- 
versible, non-equilibrium dynamics as it can appear in dissipative fluid configurations or in 
the vicinity of analogue horizons. 

In the present work we have emphasized the importance of the nut charge in the holo- 
graphic description of vorticity. In superfluids this is a quantized quantity, hence one might 
wonder whether and how nut-charge quantization could emerge in their holographic de- 
scription. The latter is in fact incomplete and the formation of vortices in rotating conden- 
sates calls for a more complete understanding, which justifies our efforts. 



24 The Cotton tensor itself is third order in the derivatives of the velocity field, but it cannot appear at this 
order in the fluid/ gravity metric because it has the wrong parity. The fourth order is indeed the smallest order 
for which it can appear. 



30 



An issue worth mentioning is the relationship of our work with alternative approaches 
of fluid/gravity correspondence. We have presented some preliminary results in Sec. |6]and 
we plan to elaborate on that subject in a forthcoming work. Many other roads seem open for 
further investigation, which we have not discussed. One could for example try to describe 
fluids in more complicated kinematic states, with multipolar vorticity - as a generalization 
of the monopole-like configurations created by nut charges, or the dipoles corresponding 
to Kerr cyclonic motions. This would require the generalization of the Weyl multipole so- 
lutions to asymptotically AdS spaces. The magnetic paradigm of the geodesic motion in 
Papapetrou-Randers geometries (see e.g. ETM ) might turn in a powerful tool for that task. 
Let us finally mention that analogue-gravity applications are very rich and diverse. Setting 
the bridge with holographic techniques would however require a more systematic study. 
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